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We characterize the set of generalized quantum measurements that can be decomposed into a 
continuous measurement process using a stream of probe qubits and a tunable interaction Hamilto¬ 
nian. Each probe in the stream interacts weakly with the target quantum system, then is measured 
projactively in a standard basis. This measurement result is used in a closed feedback loop to tune 
the interaction Hamiltonian for the next probe. The resulting evolution is a stochastic process with 
the structure of a one-dimensional random walk. To maintain this structure, and require that at 
long times the measurement outcomes be independent of the path, the allowed interaction Hamil¬ 
tonians must lie in a restricted set, such that the Hamiltonian terms on the target system form a 
Hnite dimensional Jordan algebra. This algebraic structure of the interaction Hamiltonians yields a 
large class of generalized measurements that can be continuously performed by our scheme, and we 
fully describe this set. 
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Many quantum systems either exhibit slow measure¬ 
ment read-out times or can only be probed weakly. Un¬ 
der such conditions, it is natural to monitor the systems 
continuously while simultaneously exerting some closed- 
loop feedback. Experiments can already be performed 
with such low latency that feedback can be performed 
continuously in real time [1, 2]. While generalized con¬ 
tinuous measurements have been studied [3, 4], in most 
systems the diffusive weak measurements [5] that consti¬ 
tute the continuous process must be applied via coupling 
to a probe system. Previously, we’ve studied a system 
with closed-loop feedback applied to a stream of probe 
qubits interacting with the system by a fixed Hamilto¬ 
nian [6]. Here, we investigate the possibilities that arise 
from closed-loop feedback when the interaction Hamilto¬ 
nian is itself subject to control. 

A key feature of [6] was the derivation of a reversibility 
equation which was used to restrict the class of measure¬ 
ments that admitted a continuous decomposition. This 
equation is necessary again in this work to ensure that 
the final measurement at long times is independent of 
the details of the path. Although we’ll restrict our anal¬ 
ysis to qubit probes and two-outcome measurements, we 
note that general two-outcome measurements are suffi¬ 
cient building blocks for n-outcome measurements [3, 7]. 

Consider a quantum system S undergoing a stochastic 
evolution driven by two-outcome diffusive weak measure¬ 
ments. The outcome of any particular step during the 
evolution is one of two weak measurement step operators 
M±{x). These step operators are functions of a pointer 
variable x which updates with each outcome. The ex¬ 
act feedback scheme is illustrated in Figure 1, and the 
process terminates when x reaches a fixed constant ztX. 
The reversibility condition can be written 

M^{x ± S)M±{x) (X I. (1) 



FIG. 1. The system S is continuously measured. At each 
timestep, we perform a weak measurement by preparing the 
probe \a) and tuning the interaction Hamiltonian Hpsix) 
based on a pointer variable x. The system and probe interact 
for a short time S and the probe is measured in an orthogonal 
detector basis (4>*|. The measurement result from the detec¬ 
tor is used to update the pointer variable from a; to a: ± <5 and 
the procedure is repeated with the new value. 

From the above equation, two consecutive outcomes that 
step “forward” from x to a;-|-(5 then “backward” from x+5 
to X have no net effect on |?/>). However, rather than 
track the evolution of \'tf) directly, we can express the 
total action of our procedure as the total walk operator 

[x/S\ 

M{x) = lim n M+{j5), (2) 

i=i 

For negative values of x we replace M+ with M_ above. 
We identify the endpoint operators Mi oc M{X) and 
M 2 oc M{—X) with the final measurement being decom¬ 
posed by this process. We will consider a simple model of 
probe-state interaction that will generate M±{x). In [6] 
we found that weak measurements with qubit probes 
had to form a probe-basis on the qubit Hilbert space. 
In particular, we required that the probe state, the or¬ 
thogonal quantum states of the detector, and the probe 
eigenstates of the interaction Hamiltonian, have mutu¬ 
ally orthonormal representations on the Bloch sphere. 
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For this reason, we choose the interaction Hamiltonian 
to he HPS = Fp (g) i{x), the probe state to be |0), and 
the detector states to be (±|. The operator acting on the 
system S' is e and is defined to be an a;-dependent linear 
combination of d constant Hamiltonian terms, 

d 

e{x) = ^pi{x)Hi. (3) 

The weak measurement step operators of Figure 1 then 
become 


M±{x) = 


( 4 ) 


—j=I =F -(2^)- 

V 2 V 2 2^/2 


The reversibility condition of Eq. 1 can now be rewritten 
in terms of £{x). Note that this condition need only be 
satisfied up to 0{6^) since the random walk induced on 
the pointer variable x will take 0{N‘^) steps to converge 
when N = [X/JJ. Collecting terms by orders of 5 yields 

M^{x ± 6)M±{x) = ^ + Y {dxsix) - 2£^{x)) + O(S^). 


Let a(x) be the proportionality constant in Eq. (1). We 
find that the reversibility equation reduces to 


dx£(x) = 2 £^(x) + a(x)I. (5) 


In the derivations that follow, we will ignore the a(x)I 
term, as it will not change the class of measurements that 
satisfy the reversibility equation. (In practice, the term 
can be reintroduced to help find bounded solutions.) 

Consider the set of controls that appear in Eq. (3). 
Without loss of generality, we can always assume that 
JIo = 4 since the action of I is equivalent to an overall 
phase on the probe system. The reversibility equation in 
Eq. (5) can then be rewritten as 


^d^Pk(x)ffk = ^ Pi(x)Pj(x)-. (6) 


fc =0 


ij=0 


where {•,•} is the anti-commutator. It will be useful to 
introduce the tensor FC for expressing the action of the 
anti-commutator on the matrices Hi. In particular. 


n{n-l)/2 

-{m,H,}= Y. (7) 

k=0 

We choose the matrices Hi for z > d such that they form 
a basis for Hn (C), the space of all n-dimensional complex 
Hermitian matrices. We will use F^^^ to denote the ma¬ 
trix resulting from fixing the index k. The reversibility 
equation Eq. (5) can then be read as 

/ dxPk = ^F('=)p 0<k<d . . 

\ 0 = ^F('=)p d<k. ^ ^ 


We now present our main result which characterizes 
solutions to the above equations. Let us denote F = 
spanjidi} so that e G F. We prove the following lemma 
about solutions to Eq. (8). 

Lemma 1. Any solution i(x) to Eq. (8) must lie en¬ 
tirely in V, a subspace of F that is closed under anti- 
eommutation. 

Proof. We note that if F is already closed under anti¬ 
commutation, then the reversibility equation reduces to 
an initial value problem in terms of the control coeSi- 
cients p{x) at a: = 0. However if F is not closed under 
anti-commutation, then we must characterize the set of 
vectors p such that Eq. (8) is satisfied. To do so, consider 
choosing any k > d and solving the associated equation 
pFv'^^^p = 0. Note that the matrix F*^^) is symmetric and 
defines a quadratic space over Every quadratic space 
admits a Witt decomposition [8] which in our case is 

N 

^0IF,©Ho©F'. (9) 

In the above, Wi are hyperbolic planes, Vq is the nullspace 
ofr('=), and V' is an anisotropic subspace of Solutions 
to x'^WiX = 0 are span{[l, l]}Uspan{[l, —1]}. Addition¬ 
ally, there are no vectors which satisfy SFV'x = 0 for the 
anisotropic subspace V. Let be the isomorphism 
of to and p = Then possible 

solutions to = 0 must lie in 



for a fixed choice of = ±1. To fully solve Eq. (8) we 
must now recurse the above procedure. At each step we 
restrict p to lie in the subspace V defined by a choice of 
Xi. We then define a new matrix basis for the controls 
restricted to V and generate a new set of F*^^) matrices. 
We then choose a new k and decompose V using F^^^. 
Since the order in which the k are chosen will affect the 
form of V, it is also important to enumerate all sequences 
of choices of k and Xi. This procedure terminates when 
the vector space of Hermitian matrices V formed from 
V is closed under anti-commutation. Furthermore, since 
the Witt decomposition is unique (up to isometries of 
V), we can guarantee that this procedure lists all closed 
subspaces contained in F. It remains only to show that if 
p{0) € H for a particular sequence of choices of k and Xi, 
that p{x) will remain in the same subspace for all other 
values of x. This follows directly, however, from the fact 
that if i{x) € V then e^(x) G V and so G V. □ 

Lemma 1 establishes that in order to solve the re¬ 
versibility equation, one must use a set of controls whose 
span is closed under anti-commutation. The proof of 
the lemma also includes an implicit algorithm for finding 
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closed subspaces given a set of Hermitian matrices. The 
next lemma gives the structure of the subspaces enumer¬ 
ated by lemma 1. 

Lemma 2. The e(x) operator has the form 
S(V) 

= 0 Ui{x)Di{x)Ul{x). (11) 

1=1 

where S{Y) is the number of simple components of the 
algebra V (with anti-commutation as its product), and 
Di{x) and Ui{x) correspond to the I**' simple component 
and are given by Table I. 

Proof. We begin by identifying V as a finite-dimensional 
Jordan algebra. Every such algebra accepts a 
Wedderburn-type decomposition [9, 10], 

S(V) 

( 12 ) 

1=1 

where S(V) is the number of simple components B/ of 
V. A classification of all finite-dimensional simple Jor¬ 
dan algebras was given by Jordan, von Neumann, and 
Wigner [11]. The three types of Jordan algebras that 
can be found in our decomposition are the self-adjoint 
real, complex, and quaternionic matrices. The isomor¬ 
phism in Eq. (12) leaves a lot of freedom in terms of how 
to represent each of these simple components by Hamil¬ 
tonian terms. We summarize the possible representations 
in Table I. (Note that the exceptional Albert algebra is 
absent, since octonions do not have a matrix representa¬ 
tion over M or C [12]). Since V can be written as a direct 
sum, we can also write 



TABLE I. We list all rank-n representations of Jordan alge¬ 
bras that can be embedded into a span of Hermitian matrices. 
The third representation corresponds to the 2-dimensional 
embedding of C into R. The fourth and fifth representa¬ 
tions correspond to 2- and 4-dimensional embeddings of H 
into C and R. The notation diag(R") refers to the space of 
n-dimensional diagonal matrices. 

U}{x)dxUi{x) = idxGi{x). This reduces the above equa¬ 
tion to 

i[dxGi,Di]+dxDi=2Df. (14) 

The entries of the commutator term are 

{[dxGi,Di])^. = dxgS [df^ -df^) (15) 

from which we can infer that the diagonal entries of the 
commutator term are 0 if A d[^\ We can show that 
df\x) has the solution tanh(a; — cf^), and thus for any 
i,j such that yf c'p, gf^ is constant. All together, 
Eq. (14) has the solution 

df\x) = tanh(a; — cf^) V i, 

9p\x) = gl^\0) Vi,j : c-, (16) 

9 ij{x) = gf^{x) yi,j : 


S{Y) 

£= 0ei(a;)- 

1=1 

Each operator in the direct sum can, in turn, be diago¬ 
nalized to yield the form in the statement of the lemma. 

□ 

It remains to describe the form of the endpoints of the 
continuous process Mi, M 2 . We use the reversibility and 
propagation equations to solve for them directly in the 
following lemma. 

Lemma 3. The e{x) operator and the total walk operator 
M{x) are simultaneously diagonalizable. 

Proof. We begin by noting that Eq. 5 can be solved for 
individual blocks ii (x) which yield, 

dx{Ui{x)Di{x)U}{x)^ =2{Ui{x)Di{x)U}{x)^ . (13) 

Since Ui{x) is a unitary matrix we can write it as the 
exponent of a Hermitian matrix Gi{x) and we note that 


The solution above for d!'p{x) is found by reintroducing 
the a{x)I term to our equations (which we’ve ignored 
thus far). We note that in the cases where cf ^ , g^ 

need not be constant. This, however, does not affect the 
form of £, or of Mi, M 2 . 

We now turn our attention to the total walk operator 
given in Eq. (2) which obeys the following differential 
equation (up to a normalization factor): 

dxM{x) = —s{x)M{x). (17) 

We can write M{x) in the diagonal basis of i{x) by in¬ 
troducing the operator 
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Note that since M(0) = I then iV(0) = I and any so¬ 
lution N{x) must be diagonal. Thus the total walk op¬ 
erator and the e{x) operator are diagonal in the same 
basis. □ 


Lemmas 1, 2, and 3 combined give the full characteri¬ 
zation of Ml and M 2 operators achievable by our scheme: 

Theorem 1 (Main result). A continuous measurement 
using qubit probes and closed-loop feedback on the inter¬ 
action Hamiltonian (as in Fig. 1) can realize any mea¬ 
surement {Ml, M 2 } of the form 


S(V) 


Mi = 0[// 




' rank(Ej) 

0 




Ui, 


( 20 ) 


where M 2 = {I—M\Mi)^/‘^ is diagonal in the same basis. 
The parameters ^ are real and contained in (0,1) and 
is a projector onto 1, 2, or 4 basis states. 

Proof. Recall that the number of distinct diagonal entries 
possible in Di{x) is rank (B/). However, each distinct en¬ 
try can appear 1, 2, or 4 times depending on the par¬ 
ticular representation from Table I. Using lemma 3 we 
can plug our solution for Di{x) into Eq. (19) to find that 
the diagonal entries of N{x) are exp(fg tanh(i/ —cf ^)dy). 
The total walk operator M{x) must then be 


/rank(Bi) ^ / fn\ \ 

M(a;) = [//(x) ( 0 C/i(x). 


( 21 ) 

The endpoint operators Mi and M 2 are proportional 
to M{X) and M{—X). Their diagonal entries are A^*\ 
which after renormalization approach 0 when cf^ —>■ 00 
and 1 when —>■ — 00 . □ 


Note that in theorem 1 the eigenvalues of Mi and M 2 
are restricted to lie in the open set (0,1), not the closed 
set [0,1]. This is a consequence of the reversibility con¬ 
dition at the points x = X — S and x = X 4- 6. At 
these points, setting any eigenvalue of the total walk op¬ 
erator to 0 would be effectively a projection, which is 
an irreversible operation for the random walk. However 
we can approach arbitrarily close to any such projective 
measurement. 

To allow for direct comparisons with the scheme of [6], 
we provide the following corollary. 

Corollary 1 (Spectrum of the measurement). Given the 
ability to perform any unitary transformations directly 
before and after the continuous process of theorem 1, 
one can continuously decompose any measurement with 
rank (B;) distinct singular values. 


Proof. The endpoint measurement operators Mi, M2 in 
theorem 1 can have up to rank (B;) distinct eigen¬ 

values. We can decompose any pair of general end¬ 
point operators Mi, M 2 using their polar decompositions 
M, = Wi(M/M,)i/2. Then, we can use a procedure like 
that of Figure 1 to measure the positive Hermitian oper¬ 
ators (M^Mi)^/^ and subsequently apply Wi depending 
on the measurement result. □ 

In this work we’ve characterized the full class of con¬ 
tinuous measurements achievable using a stream of probe 
qubits and a tunable interaction Hamiltonian. Given a 
set of linearly controlled Hamiltonian terms we provide 
a method to exhaustively list all continuous decomposi¬ 
tions achievable with the control set. The class we find 
has a simple block-diagonal form, but results from a non¬ 
trivial application of the reversibility condition. Notably, 
measurements in this class have a quantifiably broader 
spectrum than in the case of a fixed interaction Hamil¬ 
tonian. 

Our work makes critical use of finite-dimensional Jor¬ 
dan algebras. This is surprising since these algebras have 
had little application elsewhere in quantum mechanics. 
Our model for continuous measurements does not include 
internal dynamics Hs for the system or the probe, nor 
does it account for environment noise. In the presence of 
Hs, successive realizations of the continuous decomposi¬ 
tion would yield inconsistent results unless Hs commutes 
with the measurement operators. 

The model presented here is still not the most gen¬ 
eral description of all continuous measurements realizable 
with a stream of probes. A completely general descrip¬ 
tion would have to consider higher-dimensional probes, 
multiple outcomes to the weak measurement steps (as 
well as the endpoint measurements), and a more general 
reversibility condition. This is the subject of ongoing 
work. If Jordan algebras reappear in that scenario, then 
they will have found renewed application in quantum me¬ 
chanics. 
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